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Specific iieat, magnetic torque, and magnetization studies of LiCoP04 olivine are presented. 
They show that an unique set of physical properties of LiCoP04 leads to the appearance of features 
characteristic of 2D Ising systems near the Neel temperature, Tn = 21.6 K, and to the appearance 
of an uncommon effect of influence of magnetic field on the magnetocrystalline anisotropy. The 
latter effect manifests itself as a first-order transition, discovered at ~ 9 K, induced by magnetic 
field of 8 T. Physical nature of this transition was explained and a model describing experimental 
dependences satisfactorily was proposed. 

PACS numbers: 75.40.Cx, 75.30.Kz, 75.80.+q, 82.47.Aa 
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INTRODUCTION 

LiCoP04 olivine, crystallizing in the Prima structure^ 
Fig.[TJ exhibits a unique set of physical properties, vifhich 
makes it attractive for both basic and applied studies. 
That means: 

(i) It shovi^s an exceptionally large linear magneto- 
electric effect^""' and a large Li-ionic conductivity 
(making it promising for application as cathodes in 
Li-ions batteries))^"— 

(ii) In its structure, (100) oriented, "corrugated" Co- 
O layers can be distinguished, within which the 
Co^"*" magnetic moments are strongly coupled by 
superexchange Co-O-Co interactions. The neigh- 
boring (100) layers are coupled weakly by higher 
order interactions;^ e.g., Co-0-P-O-Co. Be- 
low T/v = 21.6 K, an antiferromagnetic ordering ap- 
pears in the system. Due to large anisotropy,®'"'^^ Co 
magnetic moments are confined to the directions ly- 
ing within the b-c plane, ca. 4.6° away from the 
b axis. The magnetoelectric effect studies^ (reveal- 
ing "butterfly" hysteresis loops and a possibility to 
produce the single domain state by application of 
a magnetic field alone), as well as the direct mag- 
netization measurementsiiii^ showed that the Co 
magnetic moments do not compensate each other 
completely and a small net magnetic moment, par- 
allel to the b axis, is present. Thus, LiCoP04 is 
an intriguing quasi - two-dimensional weakly ferro- 
magnetic Ising system. 

(iii) The presence of a spontaneous magnetization is not 
consistent with the Pnma' symmetry, for years as- 
sumed to be the magnetic symmetry of LiCoP04, 
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FIG. 1. (Color online). Orthorhombic (Pnma) olivine struc- 
ture of LiCoP04. Three unit cells (a= 10.20 A, 6 = 5.92 A, 
c = 4.70 A) stacked along the 6 axis are presented. Starting 
from the lowest one, they show, respectively, oxygen coordi- 
nations of Co (octahedral), Li (octahedral) and P (tetrahe- 
dral) ions, the examples of strong Co-O-Co and weak Co-0- 
P-O-Co superexchange couplings, and the antiferromagnetic 
ordering of magnetic moments of the Co ions. 



but is consistent with the P12]^l monoclinic sym- 
metry, in which, additionally, a nonzero dielectric 
polarization and a nonzero toroidal moment are 
allowed (the monoclinic b axis coincides with the 
pseudo-orthorhombic a axis). Attempts at measur- 
ing the dielectric polarization were unsuccessful j^ 
A nonzero toroidal moment was derive d^^'^^ on 
the microscopic level, based on the magnetic struc- 
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FIG. 2. (Color online). Specific heat of LiCoP04. (a) Temperature dependence of the total specific heat in zero magnetic 
field, B. Inset shows the A-anomaly near the Neel temperature, Tn- The parameters a = 1.31 x 10~* J/(mole K*) and b 
= —1.06 X 10~* J/ (mole K^) determine lattice contribution, (b) Magnetic contribution to the specific heat. Cm, as a function 
of temperature, measured on heating, in B parallel to the b axis. Curves for different B values are shifted along the Cm axis 
by the values given in parentheses. Inset shows the dependence of Tn on B (experimental points and fitted parabola). The 



solid lines present logarithmic dependences fitted to the experimental data near Tjv. (c) Cm — F 



ln(r),r= \T^Tn\/Tn 



Solid lines are linear approximations valid for 
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ture datai^ On the macroscopic level, four domain 
states were observed;^ two of which were inter- 
preted as " antiferromagnetic" and two other ones 
as " ferrotoroidic" . However, detailed symmetry 
considerationsi^ showed that all four domain states 
are equivalent and differ in orientation of the net 
magnetic moment. Each of the domains bears a net 
magnetic and a toroidal moment, whose signs and 
directions are mutually rigidly coupled. 

(iv) The studies of birefringence induced by magnetic 
field^'' suggest that the magnetic structure can be 
even more complex. In addition to the large, uni- 
form in space, and parallel to the b axis compo- 
nent of the main antifcromagnetic vector, L2 = 
mi — m2 — m3 -I- m4 (nii are magnetic moments 
of Co ions), small, modulated in space, perpendic- 
ular to the b axis components of L2 and of other 
antiferromagnetic vectors, defined in Ref. [Tj, can 
exist. 
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magnetoelectric,-'?'^^ transport. 
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and opticali^ii^ii^r^ properties, actual magnetic and 
electric structures of LiCoP04 and their transforma- 
tions in magnetic field have not yet been elucidated 
satisfactorily. 

Since specific heat is very sensitive to all phase transi- 
tions, this work was aimed at studying thermal properties 
of LiCoP04, at determining the order of observed phase 



transitions (spontaneous and induced by magnetic field, 
B, applied along the b axis) and at investigating how the 
intermediate dimensionality of the magnetic structure of 
LiCoP04 influences the critical behavior near T^. 



EXPERIMENT 

For the present studies, a LiCoP04 single crystal ob- 
tained by high temperature solution growth usin g lit hium 
chloride (LiCl) as flux was chosen. In Ref. [22, this 
method of crystal growth was shown to be applicable 
for the entire crystal family LiMP04 (M= Ni, Co, Fe, 
Mn). In the present case, the synthesis was realized in 
full analogy to that described in detail for LiNiP04 in 
Ref. [2J, i.e., using a molar ratio 1:3 between LiCoP04 
and LiCl in the starting mixture and using sealed plat- 
inum crucibles with 30 ml volume, with a 50 /im hole in 
the lid for equilibrating the pressure and minimizing loss 
of the highly volatile LiCl solvent. The growth parame- 
ters and the special technique for separating the fiux from 
the crystals were identical with those used for LiNiP04i^ 
No impurity phases occur in the described synthesis pro- 
cess. The growth morphology of the LiCoP04 crystals 
has been described in detail in Ref. and is character- 
ized by the development of orthorhombic (100), (210), 
(Oil), and (101) facets, which may be used as reference 
for the preparation of samples (even without X-ray ori- 
entation). 
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FIG. 3. (Color online). Magnetic contribution to tiie specific 
heat of LiCoP04, Cm, measured on heating and on cooling. 
Curves for different magnetic field values are shifted along 
the Cm axis by the values given in parentheses, (a) Lack of 
hysteresis around Tjv. (b) Thermal hysteresis near T — 9 K 
for B = 9 T. 



The specific heat of the LiCoP04 single crystal was 
measured by means of the relaxation method, using the 
Physical Property Measurement System, PPMS, made 
by Quantum Design. Estimated uncertainty of the deter- 
mined specific heat values was ~ 2% . In B = 0, studies 
were done from 2 to 300 K. Since no phase transitions ap- 
peared above Tjv, temperature dependences for nonzero 
B values, ranging from 1 up to 9 T, were measured up 
to 40 K only (the magnetic field, B, was applied along 
the b axis). The experimental points were measured ev- 
ery 0.3 K (for B^O) or 0.2 K (for B^O) below 15 K 
and every 0.1 K within the critical region around T^. In 
all figures, not all experimental points are marked with 
symbols to keep legibility. Supplementary magnetization 
and magnetic torque measurements have been performed 
by using respective measurement options of PPMS. 



RESULTS 

The zero-field temperature dependence of specific heat 
is plotted in Fig. [5^1. The inset shows the A-shaped 
anomaly accompanying the paramagnetic - weakly fer- 



TABLE I. Parameters fitting the logarithmic dependence, Eq. 
([21), to the experimental data the best. A and F are given in 
J/(mole K). For each parameter, an estimated uncertainty of 
the last digit (or of the two last digits) is given in parentheses, 
e.g., 21.63(3) means 21.63 ± 0.03. 
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romagnetic phase transition at T/y. Due to low electric 
conductivity,^'^ ~^ 10~® Scm~^, the electronic contribu- 
tion to the specific heat is negligible and the total specific 
heat, Cp, consists of the lattice, Cph, and magnetic. Cm, 
contributions only. Below 60 K, Cph can be described by 
the formula: 



Cph (T) = aT^ + hT^ 
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where the term ^ T^ represents the low-temperature de- 
pendence in the Debye model and the term ~ T^ is the 
correction of that modeP"* for the nonlinear phonon dis- 
persion relation: lo — Ci\k\ 4- C2|fcp. Based on Eq. ([1]), 
the Debye temperature was estimated to be Od = 470 K. 
The magnetic contributions to the specific heat, de- 
termined by subtracting the Cph (T) calculated accord- 
ing to Eq. H]) from the measured specific heat, as well 
as evolution of T/v with B are presented in Fig. [Jb- 
To determine the order of the phase transition occur- 
ring at T/v, Cp (T) for B ~0 and 9 T was measured on 
heating and on cooling. For both field values, no ther- 
mal hysteresis was detected. Fig. [3^, which strongly sug- 
gests that this is a second order transition and that its 
order does not change in magnetic field. It was veri- 
fied that the experimental data can not be described 
appropriately by assuming the classical form?- of criti- 
cal behavior: Cm ^ {\T ~ T]^\/Tn)^"', where the criti- 
cal exponent a takes a value between (corresponding 
to the logarithmic divergence for the two-dimensional, 
2D, Ising systero^^) and ~ 0.119 (found for the three- 
dimensional Ising modeP^). Thus, the critical behaviour 
in the form:^^ 
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was assumed and the satisfactory description of the ex- 
perimental data has been achieved, Figs.[2}D and [5};. On 
the contrary to the ideal 2D Ising system)^ the anomaly 
at T/v is evidently asymmetric with respect to T/v, Fig. 
[5^, which suggests that both A and F parameters have 
different values for both sides of Tjy. This qualitative 
expectation was confirmed by calculations and the A , 
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FIG. 4. (Color online). Field- induced first-order phase transition. Curves in panels (a) and (b) are shifted along the y-axes 
by the values given in parentheses, (a) Magnetic specific heat vs. temperature. For B ^ 8 T, the anomaly at ~ 9 K is visible. 
Vertical blue and red arrows indicate the temperatures, at which torque was measured. Solid lines present the calculated 
magnon contributions. The insets show experimental (circles) and fitted theoretical (solid lines) dependences of two parameters 
appearing in Eq. ((3]), ai (in J/(mole K^")) and 02 (in K), on B. (b) Magnetic torque for B rotating within the b-c plane {6 
is counted from the b axis). Curves plotted with full and open symbols were measured for B rotating in opposite directions. 
Dependences calculated within the proposed model (solid lines) are superimposed on the experimental curves, (c) Outline of 
the proposed two-sublattice model of the magnetic structure, (d) Superimposed magnetic torque dependences measured at 
T = 9 and 12 K. The oblique arrows indicate directions along which the torque maximum moves with increasing _B at T = 9 
K and at T = 12 K. 



F='=, and T^ values fitting the experimental data the best 
are given in Table H] (It should be stressed that the un- 
certainty of T/v given in Table U is the uncertainty of 
the theoretical, fitted parameter. For temperatures ^ 30 
K, the uncertainty of the absolute temperature values 
determined in PPMS is ±1%, whereas relative temper- 
ature changes ~ 0.03% can be detected an stabifized). 
We attribute the affect of asymmetry of the A-anomaly 
to the quasi-2D character of the magnetic structure, i.e., 
to the fact that the buckled (100) layers of strongly 
coupled Co^+ magnetic moments are not isolated but 



weakly coupled mutually. With increasing B, T/v de- 
creases parabolically, as illustrates the inset to Fig. [5]d, 
and the A-anomaly decreases, but the transition remains 
sharp. Such a behaviour is characteristic of a 2D anti- 
ferromagnetic Ising system. ^^ 

An additional anomaly appears in S = 8 T at 8.8 
K, Figs. Wfl, [3b, and Hi. For B = 9 T, it becomes more 
pronounced and shifts to 9.2 K. Near the anomaly, a hys- 
teresis between the curves measured on heating and on 
cooling appears and the anomaly measured on cooling 
is smaller. Fig. [Hb- Since these two effects are the basic 



characteristics of first-order transition o^^'^" (the first one 
is related to overheating and overcoohng phenomena and 
the second one is inherent in the relaxation method of 
measurement), we interpret the anomaly as the indica- 
tion of occurrence of a first order phase transition. It can 
be supposed that the tendency observed on increasing the 
magnetic field from 8 to 9 T will be preserved and the 
anomaly will increase and shift to higher temperatures 
with further increase of the magnetic field. 

To estimate a change of magnetic entropy related to 
this transition, we assumed that the magnon contri- 
bution, Cnia, being the only constituent of Cm apart 
from the transition, can be described in frames of the 
model developed for anisotropic antiferromagnetsi^ For 
the cases of low and high temperatures, that model pre- 
dicts, respectively: 
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where ai = ao{Ba - BY, 02 = &o(-Ba - B), oq, &o, 
and 03 are constants, Ba is a parameter of the order 
of anisotropy and exchange fields, ^b is Bohr magne- 
ton, and fcs is Boltzmann's constant. A good descrip- 
tion was achieved up to 14 K, Fig. |4^. For i? < 5 T, 
the experimental dependences can be fitted with Eq. ([3]) , 
whereas for larger fields a crossover between the behav- 
iors given by Eqs. ^ and ^ occurs at ~ 9.5 K. To 
get a satisfactory description above the crossover tem- 
perature for B = 7, 8, and 9 T, it is necessary to add 
to Eq. ((3]) a constant term, respectively, of 0.31, 0.41, 
and 0.53 J/(mole K). Tentatively, this can be ascribed 
to the fact that the model of noninteracting magnonS)^ 
within which Eqs. ([3]) and (JH) were derived, is inadequate 
close to T/v (which falls down from ~ 22 K for i? = to 
~ 18 K for J5 = 9 T). A steep changes of the calculated 
magnon contributions appearing at the transition point 
ioY B ^ 8 T, Fig. [4^, suggest that the transition is re- 
lated to a change in the stiffness of the magnon system, 
i.e. to a change of anisotropy or exchange interactions. 
After subtracting the calculated magnon contributions 
from Cm , Fig. |4^, the entropy change associated with the 
first-order transition, AS*, was calculated using the for- 
mula: AS* — Jj,^ [{Cm — C'ma) /T] dT (the temperatures 
Ti and T2 must be chosen sufficiently far from the tran- 
sition temperature, i.e., at points at which Cm = C'ma)- 
It was found AS/{kBNA) — 0.010, which is a very small 
value (e.g., the entropy change related to disappearance 
of a long-range order in a system of 1/2 spins is equal to 
ln(2) sa 0.693). 

Based on Ref. [13, where the magnetoelectric effect was 
shown to vanish in B ~ 12 T, and on Ref. |32, where de- 
struction of the antiferromagnetic ordering in the mag- 
netic field parallel to the b axis was shown to occur via 



FIG. 5. (Color online). Magnetization of LiCoP04 along the 
6 axis, measured in the magnetic field applied along the 6 axis, 
(a) Magnetization in B = 9 T (and its derivative with respect 
to T) as a function of temperature, (b) Magnetizations at 
T — 9 and 10 K (and their derivatives with respect to B) as 
a function of magnetic field. Solid line superimposed on the 
experimental dependence for 9 K was calculated within the 
two-sublattice model proposed. 



a series of spin-flip transitions, starting from i? ^ 12 T, 
we claim that the magnetic field of 9 T: (i) can influence 
the distribution of electric charges and, hence, the mag- 
netocrystalline anisotropy of LiCoP04, (ii) is too weak 
to induce a spin-flip transition. In order to verify these 
claims and to elucidate the physical nature of the discov- 
ered transition, supplementary magnetization and mag- 
netic torque measurements have been performed. 

The magnetization studies. Fig. [SI conflrmed that this 
is not a spin-flip transition, because neither in the tem- 
perature and field dependences of magnetization nor in 
their first derivatives any anomaly occurs at the transi- 
tion point. 

The magnetic torque has been measured at T = 9 and 
12 K, for the magnetic field of different value, rotating 
within the b-c plane, Figs.[4)D and[4ji. The angle 6 deter- 
mining the orientation of the field B within the b-c plane 
was counted from the b axis. The torque measurements 
showed unequivocally. Fig. |3Ji, that the transition is re- 
lated to a change in the magnetic anisotropy. This follows 
from the fact that at T = 9 K, i.e., below the transition 
point, the positions of the torque maxima evolve mono- 
tonically with increasing B (up to the highest B value of 
9 T) along the direction indicated in Fig. \^ by a blue 
oblique arrow, while at T = 12 K, i.e., above the transi- 
tion point, we observe a steep, qualitative change of the 
torque behavior. That means at T = 12 K, for i? = 6 T, 
the maxima fit into the tendency observed at T = 9 K, 
whereas for _B = 9 T, the maxima shift to the opposite 
direction than at T = 9 K. This is illustrated in Fig. HJi 
by a red oblique arrow. 



TABLE II. Parameters, which fit the best the measured dependences of torque on 9 and of magnetization along the 6 axis 
on B. Numbers in parentheses give uncertainty of the last digit (or of the two last digits) of the parameters, e.g., 21.72(4) 
means 21.72 ± 0.04. Parameters that changed as the result of the phase transition are given in bold. To is a small constant 
value (approximately, two orders of magnitude smaller than the maximum torque values measured) that must be added to the 
theoretical dependence of torque on 6, (|A.13[l . to offset a background, inherent in the measurement technique applied. 
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To describe this effect theoretically, we propose a sim- 
phfied model, Fig. |4}:, details of which are presented in 
the appendix, based on the following assumptions: 

(i) As it was suggested in Ref. 12, modulated in space, 
perpendicular to the & axis, nonzero components of 
magnetic moments exist in the magnetic structure 
of LiCoP04. 

(ii) The two-sublattice model, in which Co2 and Co3 
ions form the one sublattice, denoted as 1, and Col 
and Co4 ions form the other one, denoted as 2, can 
be apphed. 

(iii) The modulated, perpendicular to the b axis com- 
ponents average out to zero and only the net mag- 
netizations of both sublattices are essential, Fig. 

St. 

(iv) The deflection of the net magnetizations of both 
sublattices by 4.6° from the b axis and the pos- 
sibility of existence of domains, in which the sign 
of this deflection is different, can be neglected and 
it can be assumed that the net magnetizations are 
directed along the b axis. 

(iv) The sublattice magnetized "along" the field (de- 
noted as Ml in Fig. St) is rigid, has a well de- 
fined magnetization modulus, and its magnetic 
anisotropy can be analyzed by using the anisotropy 
constant Ki. 

(v) The sublattice magnetized "against" the external 
field (denoted as M2) is "weak" and behaves as an 
anisotropic paramagnet located within an effective 
field, Beff, composed of the external field, B, and 
the exchange (i.e. molecular) field, B™, produced 
by the "rigid" sublattice. 

(vi) The magnetization of the "weak" sublattice is 

equal to: M2a- = Xi<yBcSa- + X2<Tsign(Bcffa-)Scff o- + 
Xso-B^Sa^, where a identifies the components along 
the b and c axes. 

(vu) For -90° < 6* < 90°, the sublattice 1 is the rigid 
one and the sublattice 2 is the susceptible one. For 



90° < 9 < 270°, both sublattices exchange their 
behavior, that means the sublattice 1 becomes the 
susceptible one, while the sublattice 2 becomes the 
rigid one. 



(viii) For 6 — ±90° both sublattices are indistinguish- 
able, therefore the torque goes through zero at 
these angles. 

By taking the Bm, Xib, X2b, and X3b values, which fit 
the best the dependences of net magnetization on B ap- 
plied along the b axis (measured in Ref. 12 and in the 
present work. Fig. \Sjp), and by treating Ki, xic, X2c, 
and X3c as fitted parameters, a good agreement between 
the measured and the calculated angle dependences of 
torque was achieved (Fig. ^). The values of the pa- 
rameters fitting the experimental data the best are given 
in Table |lTl These values imply that at the transition 
point the sublattice magnetized "along" the field becomes 
"harder" {Ki grows), whereas the other sublattice be- 
comes "weaker" (its total susceptibility along the b axis 
remains unchanged, whereas the susceptibility along the 
c axis grows). Nevertheless, the b axis remains the easy 
magnetization direction (the sign of Ki, characterizing 
the anisotropy of the sublattice 1, does not change and 
the total susceptibility of the sublattice 2 remains larger 
along the b axis than along the c axis) . This fact explains 
why no anomalies are observed at the transition point 
on temperature and field dependences of magnetization 
measured in the magnetic field applied along the b axis. 
Fig. [5] The uncertainty of each parameter given in Table 
mi was estimated by keeping all other parameters fixed 
and checking that no noticeable change of the theoretical 
curve appears for the values of the examined parameter 
lying within the uncertainty range, whereas for the values 
beyond that range, fit quality deteriorates evidently. It 
should be mentioned that the proposed method of anal- 
ysis, in which 9 fitted parameters are involved, should be 
treated rather as a qualitative method of elucidating the 
physical processes occurring in the sample, not as an ac- 
curate method for determining physical parameters, e.g. 



CONCLUSIONS 

In conclusion, it was shown that in LiCoP04, the sec- 
ond order phase transition from the paramagnetic to 
the weakly ferromagnetic phase is accompanied by a A- 
shaped anomaly of specific heat, which can be described 
as the logarithmic divergence ([2]), characteristic of a 2D 
Ising system. The deviation from the purely 2D be- 
haviour was ascribed to the quasi-2D character of the 
magnetic structure. The first-order phase transition in- 
duced by an external magnetic field B ^ 8 T parallel 
to the b axis, appearing at ~9 K, was discovered and 
shown to be related to the change of magnetocrystalline 
anisotropy. 
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Appendix 

We assumed that the simplest description of the 
LiCoP04 antiferromagnet can be based on the two- 
sublattice approximation, in which Co2 and Co3 ions 
form the one sublattice, denoted as 1, and Col and Co4 
ions form the other one, denoted as 2 (Figs. [l]and|4l:). 
Then, we can apply a molecular field approximation, 
which must be modified in such a way that the presence 
of an extremely weak, net spontaneous magnetization, 
M°p, found experimentallyfili" will be mimicked. An 
approach to the latter effect can be based on considera- 
tions concerning a much simpler case of a ferromagnet. It 
is known that within the molecular field approximation, 
the magnetization of a ferromagnet can be determined 
graphically as the intersection point of the Brillouin func- 
tion: 



B=0 B=0.2B 
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1 , f2S 
■ coth 



1 



25* 



1 n.^ 1 
coth — y 

2S \2S" 



(A.l) 
and of the linear relation between the y parameter and 
the magnetization: 

M ksT 1 ,. , 
= — - — y B, (A.2) 

Ms g^lBSBm Bjn 

where i?m, Ms-, g, and S denote, respectively, molecu- 
lar field, saturation magnetization, (^-factor and spin of 
the magnetic ion. As shown in Fig. [SI for T = 0.5 Tq, 
in zero magnetic field the intersection point is located in 




FIG. 6. (Color online) Molecular field approximation for a 
standard ferromagnet. Red curve presents the Brillouin func- 
tion. Straight solid line presents the case of zero external field. 
The dash-dot and dashed lines show, respectively, the cases 
of the external field applied along and against the molecular 
field. 



the plateau of the Brillouin function and the spontaneous 
magnetization reaches nearly the saturation value Ms- 
Then, if the external field B ^ 0.2 Bm is applied along 
Bm, the magnetization changes only slightly, whereas if 
the same field is applied against B^m the intersection 
point shifts to the region of noticeably smaller magneti- 
zation values, where also the curvature of the Brillouin 
function and the susceptibility are larger. 

Since for LiCoP04, the temperature of the first-order 
transition is ^ 0.5 Tiq '^ 10 K, Bm (estimated based on 
the T/v value) is ^18 T, and the applied fields ranging 
from 6 to 9 T are of the order from 0.3 to 0.5 B„i, we 
can expect a similar behavior for this more complex case 
of weakly ferromagnetic antiferromagnet. 

Thus, we assumed that the sublattice 1, magnetized 
along the field, Figs.HJ; and [71 is rigid, with well defined 
modulus of its magnetization, Mq — 2toco/Kic, where 
TOco = g^J■BS = 3.26/iB is the magnetic moment of one 
Co"*"^ ion {S = 3/2) and Kic is the volume of the or- 
thorhombic unit cell (containing 4 formula units). As 
the result, magnetic anisotropy of this sublattice can be 
described by using the anisotropy constant Ki. On the 
contrary, the sublattice 2, magnetized against the field, 
is susceptible and the modulus of its magnetization is a 
function of value and direction of the field B. Thus, the 
formalism involving anisotropy constants is inapplicable 
and we can describe a magnetic anisotropy of the sublat- 
tice 2 by introducing different magnetic susceptibilities 
along different crystallographic directions. 

Additionally, the following experimental factsiiii^ 
must be taken into account: 

(i) The dependence of spontaneous magnetization, M°p, 
on temperature has a form found in Ref 12: 

M^^^{T) ^ N(0.122 - 6.5 x 10"*(r - 10.4)^), (A.3) 

where N is a coefficient needed to convert the value ex- 




FIG. 7. (Color online) Coordinate system and the symbols 
used. 



pressed in G to desired units. 

(ii) The dependence of magnetization on the field applied 
along the b axis contains terms linear and cubic in B, as 
it was found in Refflll: 



M{B) = AQ{T) + xTB + XTB^ 



(A.4) 



(Coming beyond the molecular field approximation we 
can say that the real LiCoP04 magnetic structure can 
have a form of a very weakly spread out "fan" ^ Fig. |3h. 
Then, we can imagine that the molecular field picture de- 
scribed above is a simplification of the fact that the " fan" 
of the sublattice magnetized along the filed folds slightly, 
whereas the fan of the sublattice magnetized against the 
field spreads out considerably under influence of the mag- 
netic field.) 

Besides we assumed that the deflection of the magne- 
tization by 4.6° away from the h axis and the possibility 
of existence of domains differing in sign of this deflection, 
i.e., -1-4.6° or —4.6°, has no noticeable effect on the an- 
gle dependences of torque and on the net magnetization 
value along the h axis and can be neglected. Thus, in our 
considerations was assumed that the h axis is the easy 
magnetization direction for both sublattices. Validity of 
this assumption has been verifled by direct calculations. 

Under the assumptions given above, for —90° ^ 6* ^ 
90° the magnetizations of both sublattices are given by 
the expressions: 



Afib — Mpcosai, 



Mic = Mosinofi, 



(A.5) 



M2a = XlaB^ffa + X2aSign{B^ff„)B^g„+X3aB^Sa (A.6) 

where a denotes b and c. The effective fleld, Bcfr, acting 
on the sublattice 2 consists of the applied, B, and the 
molecular, B^, fleld and is equal to: 

Befth — B cos 9 — Bm COS Qfi 

Befic = Bs\n9 — BmSvaai. (A. 7) 

Then, the free energy of the system is given by the ex- 



pression: 
F{T,B,aul 



Ki sin^ ai — AIoB cos{9 — ai) 



E ( Jxi 



cr— fc,c 



IcrB^ff^ + -X2a\Bcffa\ "f" -XScrBoffcr 



(A., 



To assure consistency with the experimental results of 
magnetization measurements, we assume that for 9 — 
and B smaller than the spin-flip field (as it is in the con- 
sidered case), also ai = and dai/dB = 0. Then, the 
theoretical resultant magnetization of the sample, calcu- 
lated according to the formula: 



M{9^0,B) 



dF 
"dB 



= Mn 



M- 



26, 



(A.9) 



should be equal to the experimental one, given by Eq. 
(|A.4|) . By comparing coefficients at different powers of B 
(in particular, the coefficient at B'^ should be equal to 0) 
we receive the following relations: 



xib ^xT + ixTBl 
X3b ^xT- 



X2b 



'ixTB 



3 ^m, 



(A.IO) 



By substituting B = into (|A.9p . we receive the equa- 
tion: 



Mo - xTBm - xTBl 






(T). 



(A.ll) 



Next, for r = 9 and 12 K, the M^^{T) parameter was 
calculated according to Eq. (JA.3P and the coefficients Xi^ 
and xf^ were determined by fitting the function (jA.4[) to 
the M{B) dependences (since we had access to the M{B) 
dependences measured for T = 9 and 10 K only, the 
Xi^ and Xa'' coefficients for T = 12 K were determined 
by extrapolating the ones found for T = 9 and 10 K). 
Then, knowing these parameters, Bm was determined by 
solving Eq. (lA.lip and the susceptibilities Xibi for i = 1, 
2, and 3, were calculated using Eq. (lA.IOp . (We found 
B„ = 20 T for T = 9 K and B^ = 17.55 T for T = 12 
K. At the first glance it seems strange that B^ varies so 
considerably, but this can be attributed to the unusual 
parabolic dependence of the spontaneous magnetization 
onTdXSl). 

In order to determine angle dependences of the mag- 
netic torque for B rotating within the b — c plane, for the 
case —90° ^ 6* ^ 90°, the ai parameter was determined 
for each 9 value by solving (numerically) the entangled 
equation: 

— — = ( I Mosinai cosai — Afo-Bsin(0 — ai) 

dai \ Mo J 

- (xifc + sign{BcSb)X2bBcSb + X36^cfffc)-SofffcBm sinai 

+ (Xlc + sign(Beffc)X2c-Bcffc + X3cB^Sc)BcScBrn COS «! 



0. 



(A.12) 



Next, Mia- and M2a values were calculated by using the 
formulae (IA.5I) - (|A.7p and the magnetic torque acting 
on the whole sample was calculated according to the for- 
mula: 

Ta = ViMx B)^ = ^^^Kc [(A/lb + M2b)B, 



(Mle + M2c)Bi, 



m Na 



K.S 



rur, 



Mosin(6' — Qfi) 



+ (xi6 + sign(Soffb)x26-Boff6 + X3b-B^ffh)Scffbsin6' 

- (xic + sign(Bcffc)x2cScffc + XzcBIsc)BcSc cos 

(A:13) 

where V, m, rrirm and Na denote, respectively, the vol- 
ume of the sample, the mass of the sample, the LiCoP04 
molar mass, and the Avogadro number. The a2 and M2 
values were determined by using the formulae: 



a2 — arctan 



fM2c 



\M, 



2b 
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Mib + Ml- (A.14) 



For the case 90° < 9 < 270°, both sublattices exchange 
their behavior, that means the sublattice 2 becomes rigid, 
whereas the sublattice 1 becomes susceptible. Thus, the 
appropriate free energy and torque values were calculated 
by using Eqs. (fO)) . (TOI) . (|XT2t . and (|XT3t . in which 
6 was replaced with 6' = 9 — 180° and the qi parameter 
determined by solving Eq. (jA.12|) was treated as a2- 

It should be mentioned that for the special values 
9 = ±90° both sublattices are indistinguishable and the 
torque is equal to zero. 

As the result of applying the procedure described 
above, in which Ki, x\ci Xic-, and X3c were treated as 
fitted parameters, the theoretical dependences of Ta on 
6*, plotted in Fig. |4)d with (black) solid lines superim- 
posed on the experimental data, have been obtained. 
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